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I. INTRODUCTION 


In classical dynamics the concept of surface of section (SOS) was introduced by Poincare 
and it has been proved to be very useful ever since [|5|. Dynamics of the Hamiltonian sys¬ 
tem with / freedoms which is a smooth diffeomorphism in 2/-dimensional phase space can 
be shown to be equivalent to discrete volume preserving mapping in (2/ — 2)-dimensional 
(sub-)phase space. Every trajectory (which is conEned to (2/ — l)-dim energy surface) 
typically crosses (2/ —2)-clim intersection of (2/ —l)-dim energy surface and (2/ —l)-dim 
SOS infinitely many times (if the system is bound and the SOS is suitably chosen). The 
so called Poincare map maps one such crossing to the next one and it is area (volume) 
preserving. Thus, for example, Hamiltonian systems of two freedoms can be reduced to 
area preserving maps of a 2-dim plane. 

There were few attempts to extend a concept of SOS reduction to quantum mechanics 
but none have been completely successful yet. The quantum Poincare map would be an 
energy-dependent propagator which would map (not yet firmly defined) amplitude dis¬ 
tribution function over the (/ — l)-dim projection of SOS to the configuration space to 
another. In bound systems, where no trajectory can escape, the quantum Poincare map 
should become unitary in the semiclassical limit (It need not be unitary in nonsemiclas- 
sical regime). Bogomolny |l| proposes semiclassical propagator for the SOS map which is 
symmetric and approximately unitary (as h —> 0). Bogomolny’s map has a clear classical 
interpretation in terms of the classical trajectories, actions along them and corresponding 
monodromy matrices, but its greatest deficiency is that it is essentially semiclassical, so 
it cannot be viewed directly as an approximation of some exact quantum propagator. 
Smilansky, Doron and Schanz |2|. [DJ propose scattering approach for diagonalization of 
bound systems. Their method is almost exact (the error is of order O (exp(—const /h)) 
because they throw away the most of atenuating (evanescent) modes with imaginary wave 
number in order to make the basis for the S'-matrix finite) but it is specialized for quantum 
billiards (Laplacians with Dirichlet boundary conditions). In fact for the billiards, Bogo- 
molny‘s method is a limiting case (as Ji —> 0) of the method of Smilansky and coworkers. 
Recently, Gutzwiller [[J] published the exact quantum surface of section map for one sys¬ 
tem, namely, the free particle on a “trombone” attached to a parallelogram. My approach 
can be viewed as a generalization of the scattering approach to arbitrary quantum Hamil¬ 
tonian system whose kinetic energy perpendicular to SOS is quadratic. Here I present a 
local version of the quantum surface of section method and prove its consistency in fairly 
rigorous manner. The central result is the proof of the decomposition formula for the 
energy-dependent Green’s function in terms of energy-dependent propagators which act 
also (or exclusively, like my quantum Poincare map) in the Hilbert space of functions over 
the configurational surface of section. Shortly after the present work had been completed 
I managed to develop a more straightforward global approach 0,1 in terms of quantum 
mechanical scattering theory, although the present local theory gives some interesting 
results which were not obtained in &0- I have also applied the present method to cal¬ 
culate the eigenenergies with the sequential number around 20 million in a nonintegrable 
system, namely the so called 2-dim semi-separable oscillator |J which dynamically has 
all the generic features even though it is geometrically somewhat special. To the best 
of my knowledge this is the only available method enabling to calculate such high-lying 
eigenstates, which demonstrates not only the conceptual importance but also the practical 
power of the present method. After this work had been completed I have been informed 
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about a related work on the exact quantum SOS method of 2—dim Hamiltonians of the 
type p 2 /2m + V(q) by Rouvinez and Smilansky |L0 . 


In the Section II we will introduce convenient notation together with all the necessary 
definitions. Then the decomposition formula for the resolvent of the energy-dependent 
Green’s function in terms of propagators which propagate from and/or onto the SOS will 
be presented and a comment on its interpretation will be given. In Section III the proof 
of the decomposition formula will be illustrated for the systems with single degree of free¬ 
dom (1-dim configuration space). The proof of the general decomposition formula will be 
given in Section IV. Additional discussion on the interpretation of the quantum surface 
of section method and conclusions are found in section V. 


II. DECOMPOSITION OF THE RESOLVENT 

First I will introduce convenient notation for posing the problem. Vectors in /-dim 
configuration space manifold C (usually C = ffl) will be denoted by q and vectors from 
the corresponding conjugate momentum space will be written as p. The discussion will be 
limited only to the case where the SOS is perpendicular to the configuration space (CS), 
SOS = {(q,p); s(q ) = 0}. More general cases can generally be transformed to the upper 
one by the use of the canonical phase space transformations. Thus, we are considering only 
configurational surface of section S 0 = {q G C\ s(q) = 0} C C. Assume that the topology 
of CS is simply connected and that function s can be globalized in a way that all members 
of the fcnnily of surfaces of section S y = {q G C\ s(q) = y} are topologically equivalent. 
Then the CS can be written as a Cartesian product C — Sx yf\, inducing the separation 
of coordinates q = ( x,y ), p = (p XJ p y ). Therefore one may write S y = (S,y). For the 
most useful example of Euclidean configuration space C with flat, mutually parallel 
surfaces of section this is already done by considering ( x , y) as coordinates in Cartesian 
coordinate system where y- axis is oriented perpendicularly to SOS. The boundary surfaces 
(x,y^) can be either at finite or infinite distance yn £ Rll {—oo, +oo}. Every surface 
of section cuts the CS in two pieces which will be referred to as upper and lower and 
denoted by the value of the index a =|, j. In arithmetic expressions, the arrows will have 
the following values, |= +1, |= —1. So, we will be considering not just one, but a whole 
family of surfaces of section at a time, where subvector x will denote the position inside 
SOS (parallel coordinates) and y will label the given SOS (perpendicular coordinate). 
This will be useful since it will turn out that the quantum SOS-SOS propagator (as will 
be dehned below) for S y uniquely determine such propagator for S y+ d y in nonlinear first 
order (Riccati-like) differential equation setting. 

My theory presented in this paper will apply to quite general class of bound Hamil¬ 
tonians (with possible generalizations to non-bound scattering problems) whose kinetic 
energy is quadratic at least perpendicularly to SOS 

H = ^P 2 y + H \ p x , x , y )- (!) 

The first term characterizes the perpendicular motion among surfaces of section and the 
second term captures the dynamics inside SOS. For some Hamiltonians which are not 
originally of this form this can be achieved by a proper choice of the SOS, i.e. the coordi¬ 
nate system (x,y). The parallel phase space coordinates (x,p x ) may not be necessarily 
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of the Weyl-type but they may belong to arbitrary Lie group manifold as will become 
more clear when the theory will be developed (see discussion in section V). 


In quantum mechanics, the observables are represented by self-adjoint operators in a 
Hilbert space hi of complex-valued functions T (q) over the CS C which obey bound¬ 
ary conditions T(;e, y-^fi = 0 and G dS or \x\ —> oo, y) = 0 and have finite norm 
fcdq\^(q)\ 2 < oo. We will use Dirac’s notation. Pure state of a physical system can be 
represented by a vector — ket | d/) which can be expanded in a convenient complete set of 
basis vectors, e.g. position eigenvectors | q) = | x,y), |\l>) = f c dq\q) (<gf| \P) = f c dq^(q)\q) 
(in a symbolic sense, since \q) are not proper vectors, but such expansions are still mean¬ 
ingful iff 'h(g) is square integrable i.e. L 2 -function). Every ket |T) G 7d has a correspond¬ 
ing vector from the dual Hilbert space 7 C, that is bra (H/1 G Li' , (T \q) = (g|T)*. Now, fix 
a given SOS, y = const. Operators x and p x can be viewed as acting on functions x ) 
of x only and therefore operating in some other, much smaller Hilbert space of square- 
integrable complex-valued functions over a given SOS S y . A family of such SOS-Hilbert 
spaces parametrized by y will be denoted by C y . Vectors in a SOS-Hilbert space C y will 
be denoted by (where subscript y will be omitted if possible). Eigenvectors of SOS 
position operators x, x\x'} = x'\x'} can provide a useful choice for the complete set of 
basis vectors of C y . {cc|^} = {fi>\x}* can be interpreted as a quantum mechanical am¬ 
plitude for a trajectory of a system being in quantum SOS-state \fii} to cross SOS S y at 
the point x. Since surfaces S y = (S,y) are all topologically equivalent, the SOS-Hilbert 
spaces C y are also all isomorphic, i.e. equivalent to some abstract SOS-Hilbert space C. 
Isomorphism I y i y G Lin(C y , C y ’) (where Lm(A,B) denotes the set of all linear mappings 
from Hilbert space A to Hilbert space B) is being provided by universal definition of the 
(symbolic) SOS-position eigenstates \x} y = \x} 


dy'y / | ® } y ' j/ {|, 


( 2 ) 


which is the well known decomposition of identity of C y for y = y'. The quantum Hamil¬ 
tonian (like the classical one ([!])) can be written as a sum of two terms 




H', H' = H'[p x ,x,y) = H'{-ihd x ,x,y ), 


( 3 ) 


where the eigenstates of the inside Hamiltonian H' restricted to the SOS-Hilbert space 
C y (treating y as a parameter and not an operator), | n} y G C y called SOS-eigenmodes 


H'\c y \n}y = E' n (y)\n} y 


( 4 ) 


provide very useful (countable n = 1,2,...) complete basis of C y as will soon become 
clear. 


Definition 1 Let us define two families of Hilbert subspaces di a EyQ C 7d, cr =|, I of L 2 - 
functions which have definite energy, i.e. they satisfy Schrodinger equation, above/below 
a given SOS S yo and are zero below/above S yo 

\ = HV{x,y) = EV(x,y) if ay > ay 0 , 

^{x,y) = 0, if ay < ay 0 }. 


4 


The values of |T) e 'E ( f yo on the SOS S yo , ^(x,y 0 ) are quite arbitrary [] and represent a 
boundary condition for the time-independent Schrodinger equation which together with 
the remaining boundary condition ^/(x,y a ) = 0 determines |T) uniquely. This suggests 
that the spaces C yo and 'H < e Vo might be of the same size (elements of both are uniquely 
determined by the values of certain function on the (/ — l)-dim SOS) and motivates 
the search for the connection between them. In order to achieve this we will use SOS¬ 
eigenmode (||) expansion of a state |T) e 'H < e vo near the surface S yo ,y = yo + ere, e > 0, 


^(x,y) = 


-im 


yo 1 


k k n (E,y 0 ) L 


out Acrki 
C n e 


■ (E,y 0 )(y-yo) _|_ c onto e -iak n {E,yo)(y-yo) 


+ 


° 0 ) • 


( 6 ) 

which is exact iff d y k n (E,y ) = 0, i.e. inside hamiltonian H 1 does not depend on the 
“parameter” y. Throughout this paper we will use the following definition of the complex 
square root \fz = + Re z) + fsgnIm ZaJ\(\z\ — Re^) which is a complex analytic 

function in upper half-plane y > 0. The accuracy of the ansatz (|6]) cannot be extended 
beyond the second order without extending the set of coefficients and differentiating quan¬ 
tities which depend on y 0 , but this is quite enough, since the present form (^) correctly 
reproduces normal derivative d y T (x, y ) and therefore also the probability currents through 
the SOS. The wavenumbers 


/ 2777 / 

k n (E, yo) — y ~jer(E — E n (y 0 )), (7) 

which correspond to SOS-eigenmodes | n} yo , assign the rest of the energy E — E n (y 0 ) 
to the perpendicular motion providing all terms in expansion (j6j) with the constant and 
correct value of the effective energy. The evanescent modes, for which E n (y 0 ) > E, have 
imaginary wavenumbers. The coefficients c° ut ,c° nto can be interpreted as the probability 
currents for SOS-eigenmodes to propagate out from the SOS and back onto the SOS, 
respectively. But only half of them are independent since the half of their nondegenerate 
linear combinations are arbitrary 


out 


I „onto 
' L 'n 


h 2 k n (E,y 0 ) 
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dx^{x, y 0 ){n\x} 


( 8 ) 


since boundary condition \I f(x,yo ) is arbitrary and the other half of their linear combina¬ 
tions are fixed 


out _ onto 


= <J\ 


-in 1 


\ mk n (E,y 0 ) 


dxd y ^>{x, y 0 ){n\x} 


( 9 ) 


since the normal derivative d y x V (x, y 0 ) is determined by the solution of the Schrodinger 
equation = E^> with boundary conditions T(c c,y 0 ) and ^{x,ya) = 0. We choose the 
coefficients c° ut to be arbitrary so the remaining coefficients c° nto are uniquely determined 
by c° ut , and make the following definition. 


1 Of course, they should decay fast enough when \x\ 


oo in case when SOS is infinite. 
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Definition 2 Take an arbitrary out-going quantum SOS-state |^ out } with the eigenmode 
expansion | , 0 out } = X]c° ut |n}. Then the propagator from SOS back to SOS (above/below 

n 

SOS if a =tl, respectively) can be defined 

T a (E, yf) e Lin(C yo , C yo ) (10) 

and the on-going quantum SOS-state |^°nt°} _ Y)cff to \n} can be written as 

n 

\r nto }^t(E,yo)\r ut }- (ii) 


Generally, 

t(E, Vo ) = < 12 ) 

n,n' 


Let us join all the information about the SOS-eigenmodes ([D and their wavenumbers in 
the definition of the waveoperator K(E,y). 

Definition 3 Linear operator K(E, y) e Lin(C y ) is by definition an operator whose spec¬ 
trum is just the set of wavenumbers {k n (E,y)} with SOS-eigenvectors \n} y 

K(E,y) = Y.K(E, y)\n) y ,{n\ = [ E ~ Tcf (13) 


With aid of this new definition and the notation of Definition one can write SOS¬ 
eigenmode expansion (|(J) more compactly ^ 




Ep [x] K-W {E ,y 0 ) 

O(|{/-j/ 0 | 2 ) 


e icrk(E,y 0 )(y-y 0 ) _|_ e ~ iry k{E ,y 0 )(y-y 0 )y 


\fj OUt } 


(14) 


Lemma 1 Expansion ( \lj\) holds for all states |T) from all spaces TL a Eyi , such that cryi < 
vyo- 

The proof follows solely from the definition [I] of the spaces H Ey , since all Tl% are equiv¬ 
alent on the domain {(ay y)\ x G S, ayo < cry < ay a }, for ayi < cry 0 . 

Definition 4 Let P rT (E,y 0 ) and Q a (E,yo) be the quantum propagators from upper/lower 
(cr =][) part of CS C to SOS S yo and vice-versa 

Qa{E, yo) e Lin(£ yo ,H ) (15) 

P*(E,y 0 ) e Lin(H,C yo ) (16) 

2 The SOS-SOS propagator T a {E,yo) can be viewed also as a scattering operator of the following 
scattering problem which is obtained from the orginal bound Hamiltonian p 2 /2m + H'(px,x,y) by 
substituting it in one part of CS, ay < ayo, by the y —flat waveguide with the Hamiltonian p 2 /2 m + 
H'(Px,x,y 0 ), for ay < ay 0 . 
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such that for arbitrary states |^°ut} ^ {^ onto | g C 


QAE,y a )\r“} 6 (17) 

{r nta \PAE,y„) e «£.„ ( 18 ) 

which are uniquely determined by the boundary condition on the SOS S yo 

(x,y 0 \Q,(E,y„)\r M } = ^Sp{x\K- 1 ' 2 (E,y 0 )(l + t,(E,y„))\r'‘ t } (19) 

{r“\P,(E,y„)\ X ,y 0 ) = k7^{V>°" to |(l + f„(E,y 0 ))K- 1/2 (E,y„)\x} (20) 


Here, perhaps, a remark is in order: Conjugated energy E* was used in fllq) in order to 
make propagator P a {E, y) complex analytic function of E rather than E*. 

Now all the necessary tools are prepared to state the main result of this paper: 


Theorem 1 The energy-dependent quantum propagator (i.e. the resolvent of the Hamil¬ 
tonian) G(E) = (E — H) 1 can be decomposed in terms of the CS-CS propagator - 
with no intersection with the SOS S y — G 0 (E,y), CS-SOS propagator P a (E,y), SOS-CS 
propagator Q a (E, y), and SOS-SOS propagator T a (E } y ) 

G(E) = G 0 (E,y) + £ Q a (E, y)(l - ± a (E, y)T a (E, y^P-^E, y) 

a 

+ ^2Q a (E,y)(l-f_ (T (E,y)f (7 (E,y))- 1 f_ a (E,y)P a ( E ,y), (21) 


and all poles of the Green’s function G(E) are in one-to-one correspondence with the singu¬ 
larities of(l-T a (E, y)T_ a (E, y )) _1 where the four propagators P a (E , y), Q a (E, y ), T a (E, y ), 
and Go(E,y) are analytic. Q 


The propagator G 0 (E,y) can be calculated from ( pl|) if the Green’s function G(E) is 
known explicitly, but it is less important, since its contribution becomes negligible when 
the energy E approaches an eigenenergy (pole). Note that (q'\Go(E, y)\q) is the proba¬ 
bility amplitude to propagate from point q to point q 1 at energy E and without hitting 
the SOS S y in between, so it is equal to zero if initial and final point, q and qf lie on 
different sides of the SOS S y . The decomposition formula fl2T|) has a very strong phys¬ 
ical interpretation. In order to see this most clearly we omit the arguments ( E , y) and 
formally expand (1 — T a T~ a )~ l in terms of the geometric series 


G = G 0 + Ql A + QfAAA + Q|T|AAA A + ... 

+ q?A + OtAAA + AAAAAA + • • • 

+ AAA + AAAAA + AAAAAAA + • • • 

+ AAA + AAAAA + AAAAAAA + • • • (22) 


The arguments (E, y) will also be omitted to increase transparency in the rest of this paper 
whenever this will cause absolutely no confusion. The meaning of this expression ( p2| ) 

3 The latter propagators are non-analytic e.g. at thresholds for opening of new modes E' n where 

k n {E' n ,y) = o. 
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when put in a sandwich between ( q'\ and | q) is that the probability amplitude to propagate 
from point q to q' at energy E (PAPPE) is a sum of: (i) PAPPE without hitting the SOS 
(Go term) plus (ii) PAPPE with single intersection with SOS (QP terms) plus (iii) PAPPE 
with double intersection with SOS (QTP terms) plus .. . The alternating symbols jj in 
a compositum like (see Figure 1) mean that the system should be on different 

sides of the SOS before and after it crosses the SOS, since typical quantum trajectory (in 
the sense of summing over paths) is continuous. The proof of the statements made in this 
section will be given in Section IV. A more lengthy discussion on the interpretation can 
be found in section V. 


III. ONE-DIMENSIONAL CASE 


First we shall prove the Theorem |T] for the special case of systems with a one degree of 
freedom / = 1. Thus the main idea of the proof will be more transparent and the reader 
will follow the next section, which is a natural continuation of the previous one, more 
easily. Ozorio de Almeida || published a similar study of one-dimensional systems in the 
context of exact quantum SOS method for separable systems. 

Now the configuration space is just the interval (or the whole real axis) C = [iji, y{\ 
and the surfaces of section S yo are single points y = yo . Thus the SOS-Hilbert spaces C yo 
are one-dimensional with one normalized basis vector 11} and the SOS-SOS propagator 
is just a simple complex-valued function T a (E,y 0 ) = {l\T a (E, y 0 )\l} which can be also 
represented with the phase shift 5 a (E,yo) between the outgoing and the incoming wave 
in the local expansion (analogous to ( [6|,[T4|) ) of the (real) wavefunction | ^E yo ) e ^ E yo 

vEy 0 {y) = ]j CQS ( k ( E > Vo)(y - Vo) + l$AE, yo)) + o (e 2 ) , (23) 

T a (E,y 0 ) = e-^(^°), (24) 

where y = y 0 +ae, e > 0. The propagators onto/from the SOS P^E, y 0 , y) = {l\P a (E, yo)\y) 
and Q a (E, y, yo) = (y\Qa(E, 2/o)|l} are now unique solutions of the Schrodinger equations 


+ k 2 (E,y))PAE,y o,y)=0 ,2, , 

(al + P(E,y))Q„(E,y,y a )= 0 


2m 


(E - V(y)) 


with boundary conditions 


(25) 


QAE,y o,»o) = P*{E,Vo,yo) = (26) 

Qo(E,y a ,yo) = P a (E,y 0 ,y a ) = 0. (27) 

It is convenient to express the SOS-CS-SOS propagators in terms of the standardized 
(real) solution of the Schrodinger equation 

QAE, yo,!/) = PAE.v.v o) = ( 28 ) 


Let the matrix elements of the Green’s functions be denoted by G(E, y', y) = {y'\G(E)\y ), 
Go(E,y 0 ,y',y) = (y'\Go(E,y 0 )\y). The decomposition formula in one dimension can now 






be written as 


G(E,y',y) 


Go (E, 7 / 0 , 7 /, y) + 

Ql(E,y',y 0 )(l - T^T^P^E.yo.y) + 
Ql(E,y',yo)(l-T l Ti)- 1 P l (E,y 0 ,y) + 
Q^E 1 y' 1 yo)(l-T l T^- 1 T l P^E,y 0 ,y) + 
Ql(E,y',y 0 )(l ~ T^T i )~ 1 T^P i (E,y 0 ,y) = 


G 0 (E,y\y)- 


-it;S 


£ (« , E W (»')i , i(!/) + ‘5 

2 sin (i( 5 ) 


(29) 

(30) 


where 5 = S^(E, y 0 ) — Si(E, yo). There are then two points to be proved: 


• Point 1 All the poles of the Green function G(E, y', y) in complex energy E plane 
should come from singularities of (1 — T^(E, yo)Ti(E, t/o)) -1 and vice versa. 


• Point 2 The residuum of the decomposed Green function G(E,y',y) at the pole 
- eigenenergy E 0 should be equal to T(t/)T( 7/ / ) where T(t/) is a normalized real 
eigenfunction at energy E 0 provided that all the propagators Go, P a , Qa and T a are 
analytic at E 0 . 


The first point is easy. The RHS of (p9| , |30| ) has a singularity if T^(E, yo)Ti(E, y 0 ) = 1. 


i.e. 


^(E,y 0 ) - Si(E,y 0 ) = 0 (mod 2 t r). 


(31) 


This happens if and only if E is an eigenvalue of the Schrodinger operator H since the 
two partial solutions of the Schrodinger equation ^% yo {y), cr =|, j (see equation (|23|)) can 
be joined into the continuous and differentiable (at SOS y = yo) eigenfunction 


dy^EyoiVo) 


dy^Ey 0 (yo) 


^ Eyo(yo) ^Eho(^o) 


(32) 


EyP 


but T^ o (t/o) = (cos5j_/cos5-f = ±1)^ (t/ 0 ). The quantization condition (|3ll,|32l) is inde¬ 
pendent of the position of the SOS y = y 0 , although in classically forbidden regions where 
the wavenumber is imaginary k 2 (E,yo) < 0 the phase-shift 8 cr (E,yo) is also imaginary 
number (plus an integer multiple of 7r) and the SOS-SOS propagator T a (E,y 0 ) becomes 
real with the magnitude which is different from unity. 

The second point is more elaborate. Let us define a global solution of the Schrodinger 
equation | v b £ ; yo ) G PL 

\^Ey 0 ) = I^Lj/o) (^ 3 ) 

where we chose the sign (± if ^5 is an even/odd multiple of n) in order to make the eigen¬ 
function ^Eoyo(y) continuous at the SOS yo- In the position representation T E yo ( x i y) = 
^ Eyo ( x , y) ± v k|r yo (cc, y), one term is always zero and the other is not. If E 0 is an eigenen¬ 
ergy which satisfies quantization condition (|3l|) then the residuum of fl30|) is 


Res G(E,y',y ) = 

Ej=EjQ 


^ E oyo (y')'& E oyo (y) 


d E S^(E 0 ,yo) — d E 5i(E 0 ,yo) 


(34) 
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In order to calculate the energy derivative of the phase shift dES a (E,y ) we should first 
study its dependence on y. Note that any wavefunction \h( y ) which satisfies the Schrodinger 
equation 

(d 2 + k 2 (E,y))A>(y) = 0, with At(y a ) = 0 (35) 

is proportional to A>^ yo (y), for ay > ayo (see lemma 0) 

| 2 > 


V(y) = a cos ( k(E,y 0 )(y - y 0 ) + \8 a (E,y Q )} + O (|y 


2/o | 


(36) 


Let us differentiate this equation with respect to y, then set y 0 = y and substitute A>{y) 
back, giving 

dy^iy) = ~k(E, y) tan (%S(E, y)) V(y) (37) 

Now, differentiate the equation ( |37| ) with respect to y, use Schrodinger equation on the 
LHS and again ( |37|) on the RHS. After cancelling the 'I' (y) on both sides one obtains 
closed first order nonlinear differential equation for the phase-shifts 


d y 5 = 2k 


d y k . 

— sm (<S), 


which is equivalent to the Riccati equation for the propagator T a 

d y T a = —2iakT c + ^(1 - T 2 ). 


The corresponding initial conditions are 


5 a (E,y a ) = 7r, T a (E,y a ) = -1 if k 2 (E,y a ) > 0, 

5 a (E,y a ) = m: - iaoo, T a (E,yJ) = 0 if k 2 (E,y a ) < 0. 

Then we define the so-called incomplete normalization constants N a (E,y ) 

y<r 


Na(E,y 0 ) = j dy 
yo 


^Eyoiv) 


-i 2 


(38) 


(39) 


(40) 


(41) 


Again we shall investigate the y —dependence of incomplete normalization constants. Since 
T E yo (y) is again proportional to some arbitrary wavefunction T (x) satisfying Schrodinger 
initial value problem (|35|) (lemma ||) we can write 


«W - j,.)) 


(42) 


Use a substitution (p42|) in the definition of N a (E,y 0 ) ([41]) , then move everything that 
does not depend on the integration variable y to the LHS and derive with respect to y 0 , 
giving 

h 2 k(E,y 0 )^ 2 (y 0 ) 


a 


yo 


-N a (E,y 0 ) 


= (Vo 


Am cos 2 (|5 ct (U,t/o)) 

After performing the differentiation and using the formulas (|37|) and (|38|) for the derivative 
of A/(y) and 5 a (E,y) with respect to y one obtains a first order differential equation for 
the incomplete normalization constants 

2 m 


dyNa + cos (8 a )N a = 

N a (E, y a ) = 0. 


(1 + cos(<f CT )), 


(43) 

(44) 
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Armed with formulas (p8| . |4d| ) we are ready to determine the energy derivative of the phase- 
shift d E 8 a (E,y) which is most valuable to us. It can be uniquely determined from the 
differential equation which expresses the uniqueness of the mixed second derivative 


d y (d E S a ) = d E (d y S a ) (45) 

since the initial conditions are also known 

d E 5 a (E,y a ) = 0. (46) 

The author has used a heuristic ansatz 


d E S a = ~N a + u(E, y) cos (5 a ) + v(E, y) sin(<5 CT ) (47) 

and found that it is consistent with the system ( fi5l , [i6| ) iff u — 0 and v = — so 

Tfi 

d E fi* = ~N a - sin(<J ff ). (48) 

Now, we can easily prove the rest of the decomposition formula with the calculation of 
the residuum ([l4|) . If E 0 is an eigenenergy then from the quantization condition m one 
has sin(5-[-(A 0 , r/ 0 )) = sin(5|(i7 0 , y 0 )) and so the denominator of ([h|) 

d E Si(E 0 ,y 0 ) — d E 5i(E 0 ,y 0 ) = —(N^(E 0 ,y 0 ) — N±(E 0 ,y 0 )) = 

2/T 

= -N{E 0 , y 0 ) = - J dy [T Eoyo (y)f (49) 

Vi 


is just a negative (complete) normalization constant which makes the residuum equal to 
the product of normalized eigenfunctions 


Re8G(E,y',y) = *(y)*{y'), *(y) 

h / =h / Q 


^Epyp ill) 

\JN(Eq, 7/0) 


(50) 


q.e.d. 


IV. PROOF OF THE GENERAL CASE 


To prove the general decomposition formula (pl|) of Theorem [T] we shall use roughly the 
same procedure as for the 1-dim case in previous section although few steps will be more 
lengthy now. 

First we have to prove (see Point 1, Section III) that all poles of the Green’s operator 
(resolvent) G(E) are due to singularities of (1 — Ti(E,y 0 )T^(E,y 0 )) in the energy region 
where the propagator T a (E, yo) is well defined. The threshold energies for opening of new 
modes E' n (y 0 ) should be excluded so that the inverse K~ X I 2 {E, y 0 ) which enters into the 
definition of T a (E,y 0 ) is well defined. 

Any wavefunction ^^(a;,?/) which satisfies Schrodinger equation above/below y 0 , 
\^E yo ) ^ H-Eyo can be represented with some |-0} e C yo in a form 


^E yo ( x ,y) = ( x ,y\Qa(E,y 0 )\^}. 


(51) 
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If Eq is such a pole, i.e. an eigenenergy of the Hamiltonian id, then the corresponding 
eigenfunction dR (x,y) and its normal derivative d y ^i(x,y) should be continuous on the 
SOS S yo . In general, the eigenenergy E 0 can sometimes be degenerate, i.e. the correspond¬ 
ing eigenspace can be more than one, say, d-dimensional spanned by the eigenfunctions 
^ n (x, y), n — 1... d. So, there should exist 2d SOS-states | an} G C yo , n — 1... d, a =|, j 


such that 


{x,y\Qj(E 0 ,y 0 )\ T n} ify>y 0 , 
(x,y\Qi(Eo,y 0 )\ | n} if y < y 0 , 


(52) 


with 

(*, yo\Qi(E 0 , y 0 )\ T n} = (x,y 0 \Qi(E 0 ,y 0 )\ j n} 
9y(x,y\Q^(E 0} y 0 )\ T n}\ y \ yo = d y (x, y\Qi(E 0 , y 0 )\ j n}\ y/yo . 


(53) 


Continuity relations ([>3|) can be rewritten by using expansion (|T4]) to calculate the normal 
derivative 


d y (x,y\Q a (E,yo)\^}\a y \a yo = cr—^{x\K 1/2 (E,y 0 ) (l -f a (E,y 0 )) |V>} (54) 

and by using the completeness of position SOS-states |a?}, 

(l + f T (S 0 ,j/o)) | T n} = (l + T^Eo^o)') | | n}, (55) 

(l -T r (£ 0 ,2/o)) I T n} = - (l ‘-f l (E 0 ,y 0 )j \ | n}. (56) 

Adding and subtracting equations (|5^j5^) one obtains important relations 

T a (E 0 ,y 0 )\crn} =\ - an}, (57) 


which mean that the operator 1 — T_ CT (i7o, y 0 )T a (E 0 , yo) is singular with |cm} being the 
corresponding right null-vectors 


(l - f-,(E 0 ,y 0 )fc(E 0 , J/o)) | crn] = 0. 


(58) 

By deriving the SOS-quantization condition ( |58D we have proved the point 1 completely, 
since the reasoning (|52[ - |58D can easily be reversed: existence of d-dim nulspace of operator 
1 — T~ a (E 0j y^Tp^Eo, y 0 ) implies the existence of d-dim eigenspace of the Hamiltonian 
H. Each d SOS-states | an},n = 1 ...d are linearly independent, since J2 n c n \an} = 0 
would imply Yln c n\^n) = 0 due to linearity of relation (|52]). The null-space of 1 — 
T_ a (E 0 , y 0 )T CT (E 0 , 7/o) which is spanned by |em}, n — 1... n is therefore also d-dimensional. 
|] We shall soon need also the complementary SOS-representation of the (complex conju¬ 
gated) eigenfunctions T*(:e,7 /) which are in analogy with (|52|) provided by propagators 
P a (E 0 , 7 /q). There exists 2d SOS-states | cm*} G C yo , n = 1... d, a = j, j, such that 


{T n*\P\(E 0 ,y 0 )\x,y) if y > y 0 , 
{in*\Pi(E 0 ,y 0 )\x,y) if y < y 0 - 


(59) 


Requiring continuity of eigenfunctions and their normal derivatives results in a sequence 
of formulas completely analogous to fl53l - |56D ending with 

{an*\f a (E 0 ,y 0 ) = {-crn*\, (60) 

4 This is another proof of the Sturm-Liouville theorem which forbids degeneracies in one dimension 
since then SOS-Hilbert space £ is 1-dimensional, and so d cannot be greater than 1. 
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(61) 


and the complementary quantization condition 

{-an*| (l - f_ a (E 0 ,yo)f a (E 0 ,yo)) = 0. 

We have still to prove (see Point 2, Section III) that the residuum of the energy de¬ 
pendent Green’s function ResE=E 0 (Q / \C'(E)\q) is equal to the sum (in case of degeneracy) 
of products of orthonormalized eigenfunctions ELi 'L n (g , )'L* (o) where complex conju¬ 
gation should be used carefully since eigenfunctions are no longer necessarily real. So let 
us assume that our sample functions which span the eigenspace at E 0 are orthonormal 
(V n \%) = Sni. The following lemma will be found very useful for calculating the residui 
of operator-valued functions. 

Lemma 2 Let A be a singular operator A G Lin(L) with d— dimensional left and right 
null-space. If \L n }, \R n } G L, n — 1... d span the left and right null space Ml = ker(yR) 
and Mr = ker(A), respectively, then 

d 

Res(s-A) _1 = J2 Inv n i[{Ln\Ri}]\Rn}{Li\ (62) 

z —U . 

n,L =1 

where Inv n i[c n i] denotes the finite-dimensional-matrix inversion, 

d d 

Inv n i[{L n \Ri}\{Li\R n A = ^{L n \Ri}Invi n '[{Li\R n ,}\ = 8 nn > (63) 

i=i i=i 

Proof Let the left and right images be denoted by X L = ML. and Xr = AC , respectively. 
We note that, 

Ir = A ft, 

n = K- 

Let us prove equation (|64| ) first. 

(i) Let |a} G Xr, so 3| b} G L such that. |a} = A\b}. Then V|c} G Ml, (c|a,} = {c|A|6} = 
0, and therefore |a} G Mr. So, Xr C Mr. 

(ii) Let |a} G Xf, so V| b} G L, {a|(A|5}) = 0, {a\A = 0, and therefore |a} G Ml So, 

X^xM l = (. Mt) L ■ 

Combining (i) and (ii) with the fact that A C B A 1 D immediately gives 

(|64[). Equation (^5j) need not be proved separately since it changes to (]64|) after formal 
substitution of the operator under consideration A —> A\ 

Let |a} G Mr = X R . So, 31im z _> 0 (2 — A)” 1 ^} = |c} since — A\c} G X R . Then, of course 
Res 2=0 (^ — A) _1 |a} = 0. Analogously for |6} G M R = X L , we have Res 2=0 {£>|(£ — A)^ 1 = 0. 
The most general form of the residuum of the resolvent which remains is G Lin(M l , Mr) , 

d 

Res(z - A) -1 = c ni\Rn}{Li\. (66) 

2=0 nf^l 

The unknown coefficients c n i can be determined by applying (^) on R, n > } and observing 
that A is null when restricted to Mr, so (z— A) _1 |i?„/} = and Res(z— A)~ l \R n '} = 

z=0 

R n i }. Since vectors |i? n /} are linearly independent, the coefficients c n i satisfy 

d 

^ ^ Cnl{Li | Rn '} &nn'i ( 67 ) 

n,l =1 


(64) 

(65) 
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so they are the elements of the inverse martrix (|63|). q.e.d. 

Let us return to physics. We have shown above (see equations ® 0) that the 
operator (1 — T^ a (E 0 , y 0 )T a (E 0 ,y 0 )) has left and right null-space of dimension d spanned 
by | — an*} and |cm}, n — 1... d, respectively. Let us now try to evaluate the residuum 
of its inverse 



Res fl -f°-(E- E 0 )d E f°) 1 = 

E=Eq \ ' 

- Res (z - (d E f°)-\ 1 - T 0 )] 1 (^T 0 )- 1 = 

z=0 \ / 



Inv n i 


{— an*\d E T°\al} 


an}{—al* 


( 68 ) 


n,l =1 


£ Invra [{T n*\d E T°\ 1 1} + {} n*\d E f°\ } l }] \an}{-al 


n,l =1 


where shorthand notation T = T_ a (E,y 0 )T a (E,y 0 ), A 0 = A(E 0j y 0 ) was introduced to 
avoid lengthy expressions. We have used lemma |2] with A = (<9£T°) _1 (1 — T°), \R n } = 

| an}, {L n | = {— an*\d E T°, differentiated T° as a product, and in the end used transfor¬ 
mations (|^7|j60| ). 

Therefore we have to study energy derivatives of the SOS-SOS propagators d E T a (E, y). 
In analogy with the procedure in Section III, we start by asking for y-dependence of 
T a (E, y) first. Take arbitrary one-side solution of the Schrodinger equation \l/(cc, y), |\R) G 
'H Eyi - By using lemma ^ we can calculate the values of the solution and its normal 
derivative on any SOS S y for ay\ < ay < ay a (i.e. everywhere < y < y^ since y\ was 
arbitrary) 


*( x,y) = Ep.{ x \K-'/*(E,y)(l+%(E,v)) m, (69) 

d y ^(x,y) = a^^{x\K 1/2 (E,y)(l~T a (E,y))\'ijj}. (70) 


Noting the completeness of position states \x} one can express | ip} from (39) and in¬ 
sert it to dZ3) and obtain very useful relation between the values and derivatives of the 
wavefunction on the SOS in terms of the SOS-SOS propagator T a 


/ 


dxd y ^>(x, y)\x} 


iaK l/2 ( 1 - f <7 )( 1 + f a y l k 1/2 j dx^(x,y)\x}. 


s 


(71) 


Differentiate the equation (|TI[ ) with respect to y and use the Schrodinger equation in a 
rather unusual form (which can be seen by the definition of the operator K(E,y )) 

(dy + K 2 (x,y)) j dx^(x,y)\x} = 0 (72) 

on the LHS, and apply formula (|7l|) again on the RHS to eliminate the normal derivatives 
d y ^(x, y). We have used d y \x} y = 0 which is compatible with the definition of the 
y —derivative of an operator-valued function which will be given in the next paragraph. 
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Since 4 >{x,y) is arbitrary for a given y, the vector / dx^>(x,y)\x} runs over entire space 
C y , so the operators in front of it on the LHS and RHS should be equal 

—at 2 = tady (at 1/2 (i - f CT )(i + f;)- 1 /? 1 / 2 )-K 1 / 2 (i-f;)(i+f CT )- 1 /?(i-f (J )(i+r CT )- 1 ^ 1/2 . 

After formal algebraic manipulation and with the definition of few new symbols 


[A,B\± = AB±BA , 


K+ = 


d y (K 1/2 ),K- 1/2 


±’ 


(73) 

(74) 


where K~ 1 / 2 = (AT 1 / 2 )" 1 , one obtains closed first order operator Riccati equation for the 
SOS-SOS propagator 


d y f a = —ia[K, T a ) + + [AA, T a ]_ + K + - T a K + T a 
with the initial condition 


(75) 


Ta(E,y„)=0 (76) 

assuming that on the boundary y = y a all SOS-wavenumbers k n (E,y a ) are imaginary. (If 
not, they can be made such if the (fictitious) hard walls on the boundaries y = y a , which 
ensure the Dirichlet boundary conditions, are infinitesimally smoothed.) 

A remark is in order concerning the derivation in the last paragraph: The derivative 
of the propagator T a (E,y) with respect to y is perfectly well defined, although T^E.y) 
and T a (E,y + dy ) act in different Hilbert spaces, since the latter are isomorphic. So for 
any such operator-valued function of y A(y) G Lin(C y ), one should define 

1 


B v A(y) = lim j (/„ +t i( y + £)/„+,,„ - U(j/)) . 


(77) 


We shall also need the //-derivatives of the SOS-CS-SOS propagators d y Q a (E,y ) and 
d y P a (E,y). Let us start by calculating the values of d y Q a (E } y)^} on the SOS S y for 
any | -0} G C y . 

(a?, y \d y Q a (E , y)\'4i}\<^y'\^y = d y ({x , y\Q a (E, y) |'0}') d y /(x } y \Q a (E, y)\'4>}\°"y'\&y = 


'im 

h 


{x\ [d y (j7- 1/2 ( 1 + T a )) - iaK^ 2 { 1 - T a )] \ijj} 


= {x,y\Q c 


-iaK - AA - K + T a 


IV>}, 


(78) 


where the equations ( [TTJ| . p74| ) and ( |75| ) and some operator algebra were applied. Note that 
the derivative of the solution of the Schrodinger equation T(£E,//) = (x, y'\Q a (E, y)\ijj} 
with respect to smooth external parameter y is again solution of the Schrodinger equation, 
d y Q a (E, y)\ijj} G H^y. The solution of the linear Schrodinger equation is a unique and 
linear function of boundary conditions, so the linear relation for boundary condition ([78]) 
holds also globally (and since ket | ijj} is arbitrary it can be omitted) 


dyQa Qo 


-iaK - AA - K+T a 


(79) 


One can completely analogously derive similar equation for another propagator P a (E,y) 


d P = 

{ - f y 1 a 


-iaK + AA - K Ii , 


(80) 
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In analogy with 1-dimensional case (although this definition slightly differs from that in 
Section III) we define also an incomplete normalization operator N a (E,y) e Lin(C y ), 

y<j 

N a {E,y) = a j dx J dy'P a (E,y)\x 1 y')(x 1 y , \Q a (E 1 y) = (81) 

5 y 

= Pcr(E, y)Q a (E, y). (82) 


We are again interested in the derivative of the incomplete normalization operator with re¬ 
spect to y. So far we have ignored the discontinuities of (x, y'\Q a (E , y) and P a (E, y)\x, y') 
at y' = y by always approaching y' —> y from the proper side, so we have to take repre¬ 
sentation (|81| ) to account for the discontinuities, and calculate 

d y N rT (E, y) = {d y P a (E,y))Q (T (E,y) + P a (E,y)d y Q (T (E,y)- 

~ cr j dxP a (E,y)\x,y)(x,y\Q C7 (E : y). 
s 


The remaining integral over the SOS can be calculated by means of initial data (|T9|,|20|) 


/ /s /s 7777 /s, a /s 

dxP a (E, y)\x, y) (x, y\Q a (E , y) = -- T (l + T a )K~\ 1 + T a ), 
s 


so after applying formulas (|79| , |80|) one can write the differential system for the incomplete 
normalization operator 


d y N a 


Na(E,y a ) 


—ia[K, N a ] + + [K _, N a }_ - f a K + N a - N a K + T (J + 

777 /s 

+ ia-^(l + T„)K-\l + T„), 

= 0, 


(83) 

(84) 


which is linear in N&, but depends again quadratically on T CT (compare with (|75|) ). Now 
we have the tools to ask for the derivative of the SOS-SOS propagator with respect to 
the energy OeT^. It can be determined as a unique solution of the first-order differential 
system 


d y (d E T„) = dE^f,,), 

a F /:,i:v j/,,1 = o, 


(85) 

( 86 ) 


which is the obvious requirement for the uniqueness of the mixed second derivative. The 
initial condition (^) is just an energy derivative of the initial condition (|76|) . The author 
has guessed the solution in analogy with 1-dimensional case (Section III) 

d E f a = -N a + ^ [k~ 2 - T a k- 2 t ) . (87) 

The reader can verify the consistency of (|87|) by inserting it into (|85|) and applying the 
formulas ([f|,|3]) and the formulas 

d E k p = 5. pk p “ 2 , 

n 

Wi- = ^\k-\k+U 
OeE = ^{[k- 2 ,k-r + d v (k- 2 )) t 
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( 88 ) 

(89) 

(90) 






which can be derived directly from the definitions and with a bit of operator algebra, of 
course. 

Now we can finish the calculation of the residuum (B3). First we use our new result 


and the formulae (57,60) and (52,59) with definition (plD to calculate the matrix 
elements 


{U*\d E T"\ V} + {ln*\d E T"\ 11} = -{Tn*|7Vf|tZ}-Un*|7Vf||Z} = 

( y\ y 

= ~S dx J 

S \y 

&nl 5 


! y \ 

dy' + J dy' \ ^* n (x,y')^i(x,y') = 


Vi 


since the eigenstates |T„) are supposed to be orthonormal ('F n |'h i ) = <5 n j. The residuum 
(|68|) is therefore very simple 


Res (1 — T) 1 = U \an}{—an* 

E=E 0 K ' i j l 

n =1 


(91) 


The calculation of the residuum of the general Green’s function (21) is now straightforward 
Res( q'\G(E)\q) = \an}{-an*\P°_ a \q} + 

<j,n 

+ T,W\Q 0 ,Wn}{-an'\^P;\ q ) = 


E=Eq 


J 2 W\Q*'W' n }W n *\P*\<i) = 

a,cr' ,n 


(92) 


This completes the proof of the general decomposition formula. 


V. DISCUSSION AND CONCLUSIONS 

In this paper I have presented and proved the decomposition of the resolvent G(E) = 
(E — E[)~ l (formula (|2T|) ) in terms of four newly defined energy-dependent propagators 
from and/or to Hilbert space of complex-valued functions over /-dim configuration space 
(CS) (kets and bras with angle brackets) to and/or from Hilbert space of complex-valued 
functions over (/ — l)-dim configurational surface of section (SOS) (kets and bras with 
curly brackets) which have the following clear physical interpretations 

• (q'\Go(E, y)\q) is a quantum probability amplitude that the system propagates from 
point q to point q' at energy E without hitting the SOS labeled by y. 

• {x\P a (E, y)\q) is a probability amplitude that the system propagates from point q 
above (if a =|) or below (if a =|) SOS to the point x on the SOS at energy E and 
without hitting the SOS in between. 

• (q\Qa{E, y)\x} is a probability amplitude that the system propagates from point x 
on the SOS to the point q in the CS above (if a =|) or below (if a =j) SOS at 
energy E without hitting the SOS in between. 
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• {x'\T a (E,y)\x} is a probability amplitude that the system propagates from point 
x on the SOS to the point x 1 on the SOS through the upper (if a =|) or lower (if 
o =|) part of CS with respect to the SOS and at energy E without hitting the SOS 
in between. 


All these propagators are analytic (holomorphic) in the upper half energy plane, Irnif > 0. 
By expressing them in terms of the resolvent of the related scattering Hamiltonian 0 one 
can show that they can have typically a finite number of poles on the real energy-axis or 
they have even no poles on the entire real energy-axis if a given SOS S y is being crossed by 
every classical trajectory, heuristically speaking (or the corresponding scattering problem, 
whoose S-matrix is just our propagator T a (E,y ) [[7|], has no bound states, rigorously 
speaking). 

Using this interpretation one can arrive at the decomposition formula (in expanded 
form (|22|)) heuristically, just by using basic postulates of quantum mechanics about addi¬ 
tion and multiplication of probability amplitudes and knowing that typical quantum path 
is continuous (but not differentiable) || (this is why the signs f and j are alternating in 
the decomposition formula (|22|).) 

In case where there is a time-reversal symmetry, the quantum Poincare (SOS-SOS) 
map is symmetric in position representation || 

{x'\f a \x} = {x\T a \x'}. (93) 



It can also be shown |10) that if the SOS-SOS operator is represented by a finite¬ 
dimensional matrix T°, n = {n'\T a \n} in a truncated basis {|n}; >0} where we discard 

all evanescent (closed) SOS-eigenmodes with imaginary wavenumbers then the matrix 
Tf, n is unitary. 

The quantum surface of section is expected to be of great practical (numerical) value 
since the dimension of the truncated matrices of the SOS-SOS propagators (which is 
O (h 1 ~^) can be orders of magnitude smaller than the dimension of the truncated matrix 

of the full Hamiltonian (which is O Thus the practical algorithms based on 

quantum SOS method for calculation of energy spectra can be developed which are orders 
of magnitude faster than the existing ones. 

The determination of the exact SOS-SOS propagator (unlike the semiclassical one as 
in |3J) is still quite a difficult task. Solving the Riccati equations (|39| , |75| ) of the propa¬ 
gator flow through the family of surfaces of section is quite elaborate, even numerically. 
One would wish to have a direct method which has been so far successful only for spe¬ 
cial systems, like free particle on a “trombone” attached to a parallelogram [Q, or Sinai 
billiard |IT|. The author succeeded to find explicit procedure for calculation of the SOS- 
SOS propagator for the so-called semi-separable systems , which are separable for y > y 0 , 
and for y < y 0 , but they are possibly discontinuous at y — y 0 . The method has already 
been applied to one such 2-dim system with generic nonintegrable classical dynamics || 
and energy levels with sequential quantum number around 20 million were easily calcu¬ 
lated (1-3 minutes of supercomputer (Convex C3860) CPU time per level). 

Some of the open problems of the subject are: (i) developing the phase-space version of 
the present theory e.g. in the Wigner-Weyl language, (ii) developing practical methods for 
evaluating exact SOS-propagators in concrete examples, (iii) establishing the connection 
between the present formalism and path integrals with geometrical constraints which are 
still to be defined for this purpose, whereas (iv) generalization of the present theory, e.g. 
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to cases of general magnetic fields where the component of the momentum perpendicular 
to SOS does no longer enter the Hamiltonian only quadratically but also linearly, and (v) 
study of the semiclassical limit were already acomplished in [|[]. 

In the present theory the algebra of the generators of parallel motions (with respect 
to the SOS) need not be just the operator algebra spanned by (x,p x ). Indeed an obvious 
generalization is possible by considering an arbitrary Lie algebra instead, since we always 
refer just to inside Hamiltonian H' (see equation (|)) as a whole and never to x and p x 
separately. For example, instead of ordinary position and momentum, x and p x , we could 
have angular momentum or spin coordinates Ji, J 2 , J 3 . 
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Figure captions 


Figure 1 Illustration of one of the terms in the expanded form of the decomposition for¬ 
mula for the energy dependent propagator (|22|) , namely Q^T^T^P^. A typical continuous 
(but nondifferentiable) path which contributes to this term is shown and the correspond¬ 
ing probability amplitudes are marked. The latter have to be multiplied and sumed 
(integrated) over the intersection points a?, x' and x" in order to yield the probability am¬ 
plitude for the entire term, namely the probability amplitude to propagate from point q 
to point q' at definite energy E and crossing the SOS exactly three times (q' \ Q \T^TiP^\q). 
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